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Abstract 

We consider a second order linear evolution equation with a dissipative term multi¬ 
plied by a time-dependent coefficient. Our aim is to design the coefficient in such a way 
that all solutions decay in time as fast as possible. 

We discover that constant coefficients do not achieve the goal, as well as time- 
dependent coefficients that are too big. On the contrary, pulsating coefficients which 
alternate big and small values in a suitable way prove to be more effective. 

Our theory applies to ordinary differential equations, systems of ordinary differential 
equations, and partial differential equations of hyperbolic type. 
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1 Introduction 

In this paper we consider abstract evolution equations of the form 

u"{t) + 2 8(t)u'(t) + Au(t) = 0, (1.1) 


with initial data 


u(0) = u 0 e D(A 1/2 ), u'(0) = «1 eH, (1.2) 

where H is a Hilbert space, A is a self-adjoint linear operator on H with dense domain 
D(A), and 5 : [0, +oo) —> [0, +oo) is a measurable function. We always assume that the 
spectrum of A is a finite set, or an unbounded increasing sequence of positive eigenvalues. 

We are interested in the decay rate of all solutions to problem (1.1)—(1.2). In par¬ 
ticular, we are interested in designing the coefficient 8(t) so that all solutions decay as 
fast as possible as t —> +oo. 

Motivation and related literature One of the motivations of the present work originates 
from control theory. Considering an equation of the form 

u"(t) + 8Bu'(t ) + Au(t) = 0, 

where B is a linear operator on H, it is known (see for example [3]) that uniform 
exponential decay of the solutions in the energy space is independent of 5 > 0. At this 
level it is already natural to seek for an optimal damping term in the class of multiples 
of B and even, more generally, of operators of the form A B where A is some symmetric 
isomorphism from H to H commuting with B. Nevertheless such optimization suffers 
from strong limitations, as we shall see below. 

In [5] an attempt was done to optimize the decay rate by a somewhat different 
method. It consists in perturbing the conservative part, namely considering an equation 
of the form 

u"(t) + 5Bu'(t ) + Au[t) + cu(t) = 0, 

where c > 0. When c and 5 tend to infinity in a certain way the decay rate can be made 
arbitrarily large. The problem here is that the equation with c large can be considered 
as driven by the operator cl rather than A: the nature of the problem is altered. Some 
of the results from [5] were later improved in [9]. 

In [2] a different strategy was used in the special case of the string equation 

Utt. 'U'xx T 8Uf 0. 

The constant dissipation 5u t was replaced by 5(x)u t , and as 8(x) approaches the 
singular potential 1/x, it was shown that the exponential decay rate can be made as 
large as prescribed. 
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In the present work we try a different approach consisting in making 8 — 8(t) time- 
dependent. This approach will turn out to be as fruitful as the ^-dependence and 
basically applicable to any self-adjoint operator A. Earlier in [4] (see also [8]), time- 
dependent damping terms of intermittent type (namely with 8(t) that vanishes on some 
intervals up to infinity) were shown to produce exponential decay provided some condi¬ 
tion on the length of the intervals of effective damping (together with the maxima and 
the minima of 8 on those intervals) is satisfied. This result was extended in [6] to some 
cases where the damping operator involves a delay term. This was one more motivation 
to examine the case of time-dependent damping. 

A simple toy model: ODEs with constant dissipation We begin our investigation by 
recalling the behavior of solutions in the simplest example where H = R and 8(t) is 
constant, so that (1.1) reduces to the ordinary differential equation 

u"(t) + 2 8u'(t) + A 2 u(t) = 0, (1.3) 

where 8 and A are positive parameters. This equation can be explicitly integrated. It 
turns out that the asymptotic behavior of solutions depends on the real part of the roots 
of the characteristic equation 

x 2 + 2 8x + A 2 = 0. (1.4) 

When 8 < A, the characteristic equation (1.4) has two complex conjugate roots with 
real part equal to —8. As a consequence, all nonzero solutions to (1.3) decay as e~ 5t . 
When 8 > A, the characteristic equation (1.4) has two real roots rq and r 2 , with 

\ 2 

rq := — 8 — V8 2 — A 2 ~ —28, r 2 := —8 + V8 2 — A 2 ~ ——. 

2o 

Every solution to (1.3) is a linear combination of e~ rit and e~ r2t . Since we are looking 
for a decay rate valid for all solutions, in this case e~ r2t is the best possible estimate, and 
it is also optimal for the generic solution to (1.3). We recall that for linear equations the 
slowest behavior is always generic, and therefore the best estimate valid for all solutions 
is sharp for the generic solution. 

A similar argument shows that for 8 = A all solutions to (1.3) are linear combinations 
of e~ xt and te~ xt , so that te~ xt is the estimate valid for all solutions in that case. 

The graph in Figure 1 represents the exponent r in the optimal decay rate e~ rt as a 
function of 8 (namely r = 8 when 8 < A and r = r 2 when 8 > A). 

It clarifies that, if we limit ourselves to constant coefficients, we cannot hope that 
all solutions to (1.3) decay better than e~ xt , and actually neither better than te~ xt . 
Moreover, beyond the threshold A, the larger is 8 the worse is the decay rate. This shows 
also that in the general setting of equation (1.1), if we restrict ourselves to constant 
damping coefficients 8{t), the best possible decay rate valid for all solutions is te~ ut , 
where u 2 is the smallest eigenvalue of the operator A. In [2] this remark was summed 
up effectively by saying that “more is not better”. 
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Figure 1: Decay rate r as a function of the constant dissipation 5 


ODEs with nonconstant dissipation The first nontrivial case we consider is the ordinary 
differential equation 

u"(t) + 2 5(t)u'(t) + A 2 u(t) = 0, (1.5) 

where now the coefficient 5(t) is time-dependent. As long as <5 e Aoc((° ; + 00 )) the 
energy of a non-trivial solution can not vanish at any finite time. Moreover, if we limit 
ourselves to coefficients 5{t) > A, once again there exists a solution which decays at 
most as te~ xt , exactly as in the case where 5(t) is constant (see Proposition 2.7). In 
other words, once again “more is not better” and the overdamping prevents a faster 
stabilization. 

Things change when we consider damping coefficients S(t ) which alternate intervals 
where they are big and intervals where they are small (below A). We obtain two results. 

• In Theorem 2.1 we prove that every exponential decay rate can be achieved through 
a periodic damping coefficient. More precisely, for every real number R there exists 
a periodic function 5(t) for which all solutions to (1.5) decay at least as e~ m . A 
possible choice for the period of 5(t) is 

i (1.6) 

hence it does not depend on R. We can also ask further requirements on S(t), 
for example being of class C°°, or taking alternatively only two values, 0 and a 
sufficiently large positive number K. 

• In Theorem 2.2 we obtain even better decay rates. Indeed, we prove that for every 
nonincreasing function p : [0, +oo) —» (0, +oo) one can design 5{t) in such a way 
that all solutions to (1.5) decay at least as p[t). In this case 5{t) is necessarily non¬ 
periodic and unbounded, but one can choose it of class C°° or piecewise constant. 
The proof of this second result relies on the first one, and the key point is that in 
the first result we can achieve any exponential decay rate e~ Rt with a coefficient 
whose period does not depend on R. 

The results for the single equation can be easily extended to systems of the form 

u'l(t) + 2S(t)u' k (t) + X 2 k u k (t) = 0 k = 1,..., n. (1.7) 
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The only difference is that in the first result (the one with a fixed exponential decay 
rate e~ m ) the period of S(t) is now 


*o = 


7 T 
2 



( 1 . 8 ) 


Once again the period is independent of R, and this is the key point to achieve 
any given decay rate through a non-periodic and unbounded coefficient. We refer to 
Theorem 2.3 and Theorem 2.4 for the details. 

The results obtained for (systems of) ordinary differential equations can be extended 
word for word to the Hilbert setting of (1.1), provided that A has a finite number of 
eigenvalues, even with infinite dimensional eigenspaces. 


PDEs with nonconstant dissipation It remains to consider the case of operators with 
an infinite number of eigenvalues. Due to our assumption on the spectrum, in this 
case H admits an orthonormal system made by eigenvectors of A, and therefore the 
evolution equation (1.1) is equivalent to a system of countably many ordinary differential 
equations. Looking at (1.8) one could naturally guess that our theory extends to the 
general setting when the series of 1/A& is convergent. This is actually true, but not 
so interesting because in most applications the series is divergent (with the notable 
exception of the beam equation, see section 5). This leads us to follow a partially 
different path. 

In Theorem 2.5 we prove once again that any exponential decay rate e~ Rt can be 
achieved through a periodic damping coefficient. The main difference is that now the 
period of S(t) is 

to = * Y 1 ( L9 ) 

A2<2(R+A!)2 k 

and hence it does depend on R. In analogy with the previous results, we can again ask 
further structure on S(t), for example being of class C°°, or taking alternatively only 
three values (instead of two). 

The fact that the period depends on R complicates the search for better decay rates. 
Our best result is stated in Theorem 2.6, where we prove that there exists 5(t) such that 
all solutions to (1.1) decay at least as a nonincreasing function p : [0, +oo) —>■ (0,1) that 
tends to zero faster than all exponentials. This universal decay rate is independent of 
the solution, but it does depend on the operator A, more precisely on its spectrum. 


More requirements on the damping coefficient The coefficients introduced in the proofs 
of the results quoted so far alternate intervals where they are close to 0 and intervals 
where they are very large. Even the coefficients of class C°° are just smooth approxima¬ 
tions of the discrete ones. On the other hand, we already know that both large values 
and values below A are needed if we want fast decay rates (see Proposition 2.7). 
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In the last part of the paper we ask ourselves whether it is essential that the coefficient 
approaches 0 or exhibits sudden oscillations between big and small values. The answer 
to both questions is negative. In the case of the ordinary differential equation (1.5) we 
show that we can achieve any exponential decay rate e~ Rt through a periodic coefficient 
5(t) which is always greater than or equal to A — £ and has Lipschitz constant equal to 
£ (where £ is a fixed parameter). Of course the period of the coefficient now depends on 
£ and R. We refer to Theorem 2.8 for the details. 

Perspectives and open problems We consider this paper as a starting point of a research 
project. Several related questions are not addressed here but could probably deserve 
future investigations. Just to give some examples, we mention hireling the optimal decay 
rates that can be achieved through damping coefficients with reasonable restrictions, 
proving or disproving that in the infinite dimensional setting there is a bound on the 
decay rate one can achieve, extending if possible some parts of the theory to operators 
with continuum spectrum, proving or disproving that random damping coefficients are 
ineffective and never better than constant ones. 

Structure of the paper This paper is organized as follows. In section 2 we state all 
our results. In section 3 we give a rough explanation of why “pulsating is better”. In 
section 4 we provide rigorous proofs. In section 5 we present some simple applications 
to partial differential equations. 

2 Statements 

For the sake of clarity we present our results in increasing order of complexity. We start 
with ordinary differential equations, we continue with systems of ordinary differential 
equations, and finally we consider the more general Hilbert setting. In the last sub¬ 
section we investigate the same problems with additional constraints on the damping 
coefficients. 

In the sequel (t — fo) + stands for rnax{f — t 0 , 0}. 

2.1 Ordinary differential equations 

To begin with, we consider the ordinary differential equation (1.5). In the first result 
we achieve any given exponential decay through a periodic damping coefficient. 

Theorem 2.1 (Single ODE, fixed exponential decay rate). Let A and R be positive real 
numbers, and let t 0 be defined by (1.6). 

Then there exists a t 0 -periodic damping coefficient 5 : [0, +oo) —» [0, +oo) (which 
one can choose either of class C°° or piecewise constant) such that every solution u(t ) 
to (1.5) satisfies 

\u'(t)\ 2 + A 2 |w(t)| 2 < (|u'(0)| 2 + A 2 |m( 0)| 2 ) exp (— R(t — i 0 ) + ) Vi > 0. (2.1) 
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The typical profile of a damping coefficient realizing a fixed exponential decay rate 
is shown in Figure 2 on the left. In each period there are two impulses of suitable height 
K and duration p, one at the beginning and one at the end of the period. For the rest of 
the time the damping coefficient vanishes, which means that there is no dissipation. Of 
course, when the coefficient is extended by periodicity, the impulse at the end of each 
period continues with the impulse at the beginning of next period, thus giving rise to a 
single impulse with double time-length. 

These piecewise constant coefficients with only two values sound good for applica¬ 
tions. In any case, since (1.5) is stable under L 2 perturbations of the coefficient (see 
Lemma 4.1), the same effect can be achieved through a smooth approximation of the 
piecewise constant coefficient. 

In the second result we show that every fixed decay rate can be achieved if we are 
allowed to exploit non-periodic and unbounded damping coefficients. 

Theorem 2.2 (Single ODE, any given decay rate). Let X be a positive real number, let 
tp : [0, Too) —> (0, Too) be a nonincreasing function, and let to be defined by (1.6). 

Then there exists a damping coefficient 5 : [0, +oo) —>■ [0, +oo) (which one can choose 
either of class C°° or piecewise constant) such that every solution u(t ) to (1.5) satisfies 

\u\t)\ 2 + X 2 \u(t)\ 2 < (K(0)| 2 + A 2 H0)| 2 ) -<p(t) Vi > t 0 - (2.2) 

The typical profile of a damping coefficient realizing a given decay rate in Theo¬ 
rem 2.2 is shown in Figure 2 on the right. If consists in a sequence of blocks of the 
type described after Theorem 2.1, the only difference being that now the values of the 
parameters K and p are different in each block. 




Figure 2: possible profiles of 5(t) in Theorem 2.1 (left) and Theorem 2.2 (right) 


2.2 Systems of ordinary differential equations 

The results for a single ordinary differential equation can be extended to systems. The 
following statement is the generalization of Theorem 2.1. 
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Theorem 2.3 (System of ODEs, fixed exponential decay rate). Let n be a positive 
integer, let (Ai,..., A n ) G (0, Too) n , let R be a positive real number, and let t 0 be defined 
by (1.8). 

Then there exists a t^-periodic damping coefficient 5 : [0, +oo) —» [0, Too) (which 
one can choose either of class C°° or piecewise constant) such that every solution 
(ufit ),..., u n {t)) to system (1.7) satisfies 


£(WWI 2 + ^Mt)| 2 )< 

k =1 


E(W(°)l 2 + A *M0)l 2 ) 

_k =1 


exp (~R(t - t 0 ) + ) 


for every t > 0. 


The typical profile of a damping coefficient realizing a fixed exponential decay rate 
for a system is shown in Figure 3 on the left. Now the period [0,io] is the union of k 
subintervals (k — 3 in the figure), where the i-th subinterval has length 7r/(2Aj). In 
each subinterval we exploit once again the profile with two impulses described after 
Theorem 2.1. We can assume that in all subintervals the values of the parameters K 
and p are the same (and hence the time-length of the vanishing phase is different), so 
we end up once again with a damping coefficient with just two values. 

The following statement is the generalization of Theorem 2.2 to systems. 

Theorem 2.4 (System of ODEs, any given decay rate). Let n be a positive integer, let 
(Ai,...,A n ) G (0, Too) n , let ip : [0, +oo) —* (0, +oo) be a nonincreasing function, and 
let to be defined by (1.8). 

Then there exists a damping coefficient 5 : [0, +oo) —> [0, +oo) (which one can choose 
either of class C°° or piecewise constant) such that every solution (ufit),... ,u n (t )) to 
system (1.7) satisfies 


E(KMI 2 + A2Mt)| 2 )< 

k =1 


E(K(o )| 2 + aJK(o)| 2 ) 

_k =1 


■ip(t) 


Vt > t 0 . 


As in the case of a single equation, the typical profile of a damping coefficient realizing 
a given decay rate for a system is a sequence of blocks of the same type used in order 
to realize exponential decay rates for the same system, just with different values of the 
parameters K and p in different blocks (see Figure 3 on the right). 


2.3 Partial differential equations 

We examine now equation (1.1) in the general Hilbert setting. As usual, we consider 
weak solutions with regularity 

u G C° ([0, +oo), D(A 1/2 )) 0 C^QO, Too) , H). 

In the first result we achieve once again a given exponential decay rate through a 
periodic damping coefficient. In contrast with the case of ordinary differential equations 
or systems, the period of the coefficient now does depend on the decay rate. 
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Figure 3: possible profiles of S(t) in Theorem 2.3 (left) and Theorem 2.4 (right) 

Theorem 2.5 (PDE, fixed exponential decay rate). Let H be a Hilbert space, and let 
A be a self-adjoint nonnegative operator on H with dense domain D(A). Let us assume 
that the spectrum of A is an increasing unbounded sequence of positive real numbers 
{A^>r (with the agreement that A& > 0 for every k> 1). 

Let R be a positive real number, and let t 0 be defined by (1.9). 

Then there exists a to-periodic damping coefficient 5 : [0, +oo) —> [0, +oo) (which 
one can choose either of class C°° or piecewise constant) such that every weak solution 
to equation (1.1) satisfies 

\u\t)\ 2 + \A 1/2 u(t)\ 2 < (|n'(0)| 2 + \A 1/2 u(0)\ 2 ) exp (—R(t - t 0 ) + ) (2.3) 

for every t > 0. 

The typical profile of a damping coefficient realizing a fixed exponential decay rate 
for the full equation (1.1) is shown in Figure 4. Now the period [0, to] is divided into two 
subintervals of the same length to/2. In the second subinterval the damping coefficient 
is constant. The first half of the period is in turn divided into subintervals of length 
7r/(2Aj), where the Aj’s are those which contribute to the sum in the right-hand side of 
(1.9). In each of these subintervals we have again the same profile with two impulses 
described after Theorem 2.1. We can assume that the parameters K and p are the same 
in all subintervals of the first half of the period, but the constant in the second half of 
the period might differ from K. Therefore, now the damping coefficient takes in general 
three values instead of two. 

In the second result as usual we allow ourselves to exploit non-periodic and un¬ 
bounded damping coefficients. What we obtain is a decay rate which is faster than all 
exponentials. This decay rate does depend on the operator. 

Theorem 2.6 (PDE, decay rate faster than all exponentials). Let the Hilbert space H 
and the operator A be as in Theorem 2.5. 

Then there exist a damping coefficient 5 : [0, +oo) —* [0, +oo) (which one can choose 
either of class C°° or piecewise constant) and a nonincreasing function <p : [0, +oo) —> 
(0, +oo) such that 

lim ip(t)e m = 

t —I - OO 
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0 Vi? > 0, 


(2.4) 























Figure 4: possible profiles of 5(t) in Theorem 2.5 

and such that every weak solution to equation (1.1) satisfies 

\u'(t)\ 2 + \A l/2 u(t)\ 2 < O'(0)| 2 + |A 1 / 2 u(0)| 2 ) -ip(t) Vt> 0. (2.5) 

As in the case of ordinary differential equations or systems, the standard profile of a 
damping coefficient provided by Theorem 2.6 is a sequence of blocks of the same type as 
those introduced for Theorem 2.5. The main difference is that now also the time-length 
to is different in different blocks, and increases with time. The reason is that now to 
depends on R , and when R increases a larger number of eigenvalues contributes to (1.9), 
and hence to increases as well. 

2.4 Further requirements on the damping coefficient 

In the last part of the paper we investigate which features of the damping coefficient 
are essential when one wants to achieve fast decay rates. To begin with, in the follow¬ 
ing statement we list three simple situations where the decay rate of solutions can be 
bounded from below. We point out that, when looking for estimates from below, there is 
almost no loss of generality in considering just the ordinary differential equation (1.5). 

Proposition 2.7 (Estimates of the decay rate from below). Let us consider equa¬ 
tion (1.5) for some positive real number A and some nonnegative damping coefficient 

$ e ^Lc((0) +°°))- 

(1) If 5 G L 1 ((0, Too)), then all nonzero solutions do not decay to zero. 

(2) If there exists a constant M > 0 such that 5(f) < M for every t > 0, then all 
solutions satisfy 

\u\t)\ 2 + X 2 \u(t)\ 2 > O'(0)| 2 + A 2 H0)| 2 ) e~ AMt Vt > 0. 
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(3) If there exists T > 0 such that 8(t) > A for every t > T, then there exist T* > T 
and a solution u{t ) to (1.5) such that 

\u(t)\ > te~ xt Vt > T*. (2.6) 

As a consequence, if we want all solutions to (1.5) to decay faster than a given 
exponential, we are forced to choose a damping coefficient 5(t) which alternates intervals 
where it is large enough, and intervals where it is smaller than A. If we want solutions 
to decay faster than all exponentials, we also need 8(t) to be unbounded. 

In some sense these are the unique essential features. In the following result we 
show that any fixed exponential decay rate can be achieved through a periodic damping 
coefficient which is greater than or equal to A — e, and has Lipschitz constant equal to 
e, and hence it exhibits very slow transitions from small to large values. 

Theorem 2.8 (Single ODE, with borderline constraints on 5(t)). Let A and R be positive 
real numbers, and let e G (0, A). Then there exist a positive real number to, with 

t 0 <16 (- + l) R+ 2<y7T + 1 ^ + y + 2 + 8 log2, (2.7) 

and a t 0 -periodic function 5 : [0, +oo) —> [0, +oo) such that 

|<5(£) — <5(s)| < e\t — s\ Vf > 0 Vs > 0, 

5(t) > A — e Vt > 0, 
and such that every solution to (1.5) satisfies 

|-u'(f)| 2 + A 2 |n(t)| 2 < (^'(O)! 2 + A 2 |n(0)| 2 ) exp (—R(t — t 0 ) + ) Vt > 0. (2.8) 

Theorem 2.8 deals with the case of a single ordinary differential equation, but an 
analogous result holds true also for systems or the abstract equation (1.1). In order to 
contain this paper in a reasonable length we spare the reader from the details. 

More delicate is achieving decay rates faster than all exponentials through damping 
coefficients with small Lipschitz constant. This could be done at most in the same spirit 
of Theorem 2.6, the reason being as usual that now t 0 depends on R. We do not address 
this issue in this paper. 

3 Heuristics 

In this section we provide an informal description of the strategy of our proofs. Our aim 
is clarifying why pulsating damping coefficients are more effective when we are interested 
in damping all solutions to an equation or system. 

Let us start with the single ordinary differential equation (1.5). We want to design 
8(t) so that all solutions decay as fast as possible. The first naive idea is to choose 8(t) 
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very large. Bending the rules a little bit, we can even imagine to choose a damping 
coefficient which is not a function, but a Dirac delta function (which is actually a 
measure) concentrated at time t — 0, or even better a delta function multiplied by a 
large enough constant k. 

An easy calculation shows that such an extreme damping has a great effect on the 
solution with initial data tt(0) = 0 and u'(0) = 1, whose energy is instantly reduced by a 
factor e~ k . On the contrary, it has no effect on the orthogonal solution with initial data 
■u(O) = 1 and u'(0) = 0. This can be explained by observing that the damping coefficient 
multiplies u'(t), and in the case of the second solution this time-derivative vanishes when 
the delta function acts. The effect on any other solution is a linear combination of the 
two, namely highly reducing the time-derivative but leaving the function untouched. 

This apparently inconclusive approach suggests a first strategy: if we want to dampen 
a single solution, we can use a delta function acting when the energy of the solution 
is concentrated on the time-derivative. So we take again the solution with initial data 
■u(O) = 1 and u'(0) = 0, and we apply no dissipation until u(t) vanish. This happens for 
the first time when t = 7t/( 2A). At that point we apply a second delta function. 

Summing up, a first delta function at time t = 0 cuts the first solution, then the 
damping coefficient vanishes until the second delta function at time t = 7t/( 2A) cuts 
the second solution. What happens to all other solutions? Since the equation is linear, 
cutting two linearly independent solutions is equivalent to cutting all solutions. If we 
repeat this procedure by periodicity, we can achieve any exponential decay rate. If at 
each reiteration we increase the multiplicative constant in front of the delta functions, 
we can achieve any given decay rate. This is the idea behind the proofs of the results 
stated in section 2.1, and the point where the special time (1.6) and the profiles of 
Figure 2 come into play. 

The idea for systems is a simple generalization. We use a first block of two delta 
functions with time-gap of tx/{2\i) in order to dampen solutions of the first equation, 
then we use two more delta functions with time-gap of 7 t/( 2A2) in order to dampen 
solutions of the second equation, and so on. In other words, we take care of the equations 
of the system one by one. We end up with the special time (1.8) and the profiles of 
Figure 3. We are quite skeptic about the possibility of reducing the time (1.8), unless 
Aj’s satisfy special rationality conditions. 

When we deal with a partial differential equation, which we regard as a system of 
countably many ordinary differential equations, we exploit a mixed strategy. If we want 
to achieve a given exponential decay rate, a suitable constant damping coefficient does 
the job for all components corresponding to large enough eigenvalues. Thus we are left 
with cutting a finite number of components, and this can be done as in the case of 
finite systems. As a consequence, now a good damping coefficient consists in a constant 
damping half the time, alternated with a train of delta functions in the remaining half 
of the time. This is the idea behind the proof of Theorem 2.5 and the profile of Figure 4. 

In the case of partial differential equations, things are more complex if we want 
to achieve a decay rate faster than all exponentials. Indeed, when we reiterate the 
procedure, a better decay requires more components to be treated separately, and in 
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turn this implies a longer wait. The compromise between faster decay rates and longer 
waiting times gives rise to the operator dependent rate (pit) of Theorem 2.6. 

This discussion motivates also the last part of the paper. Indeed, a train of delta 
functions (or a suitable approximation) emerged as a common pattern of the damping 
coefficients which realize fast decay rates. In a first stage this led us to suspect that a 
bound on the Lipschitz constant of the coefficient, or the impossibility to attain values 
close to zero, could yield a bound from below on the decay rate of solutions. 

In Theorem 2.8 we show that this is not the case, because any exponential decay 
rate can be realized through a damping coefficient 8(t) with arbitrarily small Lipschitz 
constant. The construction of 8{t) is more involved, but once again it acts in two steps. 
In a first phase the coefficient grows and kills a first solution. Then the coefficients goes 
below A and stays there until the orthogonal solution has rotated enough so that it is 
ready to be damped by a second growth of the coefficient. 

This two-phase action (destroy the first solution, wait for rotation, destroy the second 
solution) seems to be the quintessence of all the story. 

4 Proofs 

In this section we prove our main results, following the same scheme of the statement 
section. We begin by investigating how solutions to (1.1) depend on the damping coef¬ 
ficient. 

Lemma 4.1 (Continuous dependence on the damping coefficient). Let H be a Hilbert 
space, let A be a self-adjoint nonnegative linear operator on H with dense domain D(A), 
and let T be a positive real number. Let <5i : [0, T] —> [0, +oo) and d 2 : [0, T] [0, +oo) 
be two bounded measurable functions. Let u\{t) and u 2 {t ) be the solutions to (1.1) with 
8(t) replaced by 5i(t) and 5 2 {t), respectively, and with initial data 

■Ui(O) = u 2 (0) = Uq E D(A l t 2 ), M i(0) = ^4(0) = Ui E H. 

Then for every t E [0, T] the following estimate holds true 

\u' 2 (t) - wi(t)| 2 + \A 1/2 (u 2 {t) - ui(t))| 2 < 2 (m 2 + \A 1/2 u 0 \ 2 ) e 2t ■ 

|<S 2 (s) - <5i(s)| 2 ds. (4.1) 



Proof To begin with, we observe that 

K(f)| 2 + \A^ 2 Ul (t)\ 2 < M 2 + |^ 1/2 « 0 | 2 Vt E [0, T]. (4.2) 

This inequality holds true because the left-hand side is a nonincreasing function of 
time. Now let us set 

E{t) := | u' 2 (t) - w'i(t)| 2 + | A 1/2 (u 2 {t) - Mi(t))| 2 . 
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An easy computation shows that 

E'(t) = -AS 2 (t)\u' 2 (t) -v^{t)\ 2 + A{S 1 {t)) -5 2 (t)) • «(t),t4(t) -u[(t)). 

The first term in the right-hand side is less than or equal to zero. Keeping (4.2) into 
account, we can estimate the second term and obtain that 

E'(t) < 2\u' 2 (t) — u[(t)\ 2 + 2\u((t)\ 2 • |^ 2 (t) — di(f)| 2 

< 2 Eft) + 2 (|ui| 2 + |A 1/,2 u 0 | 2 ) \$2(t) — <*i(*)h 

Integrating this differential inequality, and recalling that E( 0) = 0 because the initial 
conditions of u 2 (t) and ufit) are the same, we conclude that 

Eft) < 2 (|«i| 2 + |A 1/2 u 0 | 2 ) e 2t f |5 2 (s) - <b(s)| 2 ds Vfe[0,T], 

Jo 

which proves (4.1). □ 


4.1 Ordinary differential equations 


The following result is the fundamental tool in our theory. 

Lemma 4.2 (Decay for two orthogonal solutions to a single ODE). Let A and M be 
positive real numbers, and let t 0 be defined by (1.6). For every positive integer n, let us 
consider the function 5 n : [0, t 0 \ —* [0, +oo) defined by 

Mn if t e [0,1 /n\ U [t 0 — l/n,to], 

0 otherwise , 


5 n (t) := 


(4.3) 


and the differential equation 


u"{t ) + 2 5 n (t)u'(t) + A 2 u(t) = 0. 


(4.4) 


Let v n (t) be the solution with initial data v n (0) = 0 and i^(0) 
solution with initial data u; n (0) = 1/A and w/(0) = 0. 

Then 


limsup (\v' n (t 0 )\ 2 + A 2 |v n (t 0 )| 2 ) <e 4M , 

n —>-+oo 


limsup(|w/(t 0 )| 2 + A 2 |w„,(t 0 )| 2 ) < e 4M . 

n—>•+oo 


1. Let w n (t ) be the 


(4.5) 

(4.6) 


Proof For every solution u(t) to (4.4), let us consider its energy 

E u {t) := \u'{t)\ 2 + A 2 |u(t)| 2 . 

A simple computation of the time-derivative shows that E u (t) is nonincreasing. Let 
us set for simplicity 

t n := 1/n, s n := t 0 - 1/n, 

and let us assume that n is large enough to that t n < s n , and hence the two intervals 
where 6 n (t) = Mn are disjoint. 
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Estimate on v n (t) Since E Vn (t ) < E Vn (0) = 1, it follows that \v' n (t)\ < 1 for every 
t G [0, t 0 ]- Thus from the mean value theorem we obtain that |v n (0) — v n (t n )\ < t n , and 
hence 

lim v n (t n ) = 0. (4.7) 

n—>•+ oo 

Let us consider now the time-derivative. To this end, we interpret (4.4) as a first 
order linear equation in u'(t), with forcing term —A 2 u(t). Integrating this differential 
equation we obtain that 

v' n (t) = v' n (0)e- 2Mnt + [ e 2Mn( ' s ~ t ' > X 2 v n (s) ds Vt G [0,t n ]. 

Jo 

Now we set t = t n , we recall that v' n {tS) = 1, and we pass to the limit as n —> +cx). 
The integrand is bounded because s < t and |u n (s)| is bounded owing to the energy 
estimate. Since t n —> 0 the integral tends to 0, and hence 

lim v' n {t n ) = e~ 2M . (4.8) 

n— >-+oo 

From (4.7) and (4.8) we conclude that 

lim (\v' n (t n )\ 2 + X 2 \v n (t n )\ 2 ) =e~ iM , 

n—H-oo 

which implies (4.5) because the energy is nonincreasing with time. 


Estimate on w n (t ) Let us begin with the interval [0, t n ]. The same argument ex¬ 
ploited in the case of v n (t) now leads to 


lim w n (t n ) = 1/A, lim w' n (t n ) = 0. (4.9) 

n^-H-oo n^-+oo 


Roughly speaking, this means that nothing changes for w n (t) in the interval [0, l n ]. 
Let us consider now the interval (t n ,s n ), where 5 n (t) is identically 0. An easy com¬ 
putation shows that in this interval w n (t) is given by the explicit formula 

in f (f ^ 

w n (t) = w n (t n ) cos(A (t - t n )) + n sin (A (t - t n )). 

Setting t = s n we find that 


w n (s n ) 


w n (t n ) cos 




+ 



sin 




5 


w' n (s n ) = -Xw n (t n ) sin 


7T 2A 


n 


+ w' n (t n ) cos 


7T 2A 
2 n 


Passing to the limit as n —> +cxd, and keeping (4.9) into account, we deduce that 


lim w n (s n ) = 0, 

n—>•+oo 


lim w' n (s n ) = -1. 

n^-H-oo 
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Roughly speaking, this means that the interval [t n , s n ] has produced a rotation of 
w n {t) in the phase space, with the effect of moving all the energy on the derivative. 

Let us finally consider the interval [s n , to], where we argue as we did in [0, t n \ with the 
function v n (t). Due to the uniform bound on w' n (t) coming from the energy estimate, 
from the mean value theorem we deduce that 

lim w n (t 0 ) = lim w n (s n ) = 0. (4.10) 

n—>■+oo n—>-+oo 

As for the derivative, once again we interpret (4.4) as a first order linear equation in 
u'(t), and we find that 

w' n {t) = w' n {s n )e- 2Mn{t ~ Sn) + r e 2Mn ^X 2 w n (s) ds Vt G [s n , t 0 ]. 

J S n 

Once again the integrand is bounded because s <t and w n (s) is uniformly bounded 
owing to the energy estimate. Setting t = t 0 , and passing to the limit as n —> +oo, we 
conclude that 

lim v/ n (to) = —e~ 2M . (4.11) 

n—>•+oo 

From (4.10) and (4.11) we deduce (4.6). □ 

A careful inspection of the proof reveals that the inequality in (4.6) is actually an 
equality, and the limsup is actually a limit. With similar arguments one could show 
that the same is true in (4.5) (it is enough to follow the solution v n (t) until the end of 
the interval). In any case, (4.5) and (4.6) are what we need in the sequel. 

In the next result we apply Lemma 4.2 in order to dampen all solutions to (1.5). 

Lemma 4.3 (Decay for all solutions to a single ODE). Let A and M be positive real 
numbers, and let to be defined by (1.6). 

Then there exists a bounded measurable damping coefficient 5 : [0,f 0 ] —> [0, +oo) 
such that every solution to (1.5) satisfies 

K(t 0 )| 2 + A 2 HO)| 2 < (K(0)| 2 + A 2 H0)| 2 ) • 2e~ M . (4.12) 

Furthermore, one can choose 5{t ) such that 

• either it is of the form (4-3) for a large enough n, 

• or it is of class C°° and all its time-derivatives vanish both int = 0 and int = t 0 . 

Proof We claim that (4.12) holds true for all solutions to (1.5) if the damping coefficient 
is of the form (4.3) with n large enough. To this end, let 5 n (t), v n (t) and w n (t) be defined 
as in Lemma 4.2. Due to (4.5) and (4.6), the two inequalities 

K(t 0 )| 2 + A 2 K(0)| 2 < e~ M , \w' n (t 0 )| 2 + X 2 \w n (t 0 )\ 2 < e~ M (4.13) 
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hold true provided that n is large enough. Every solution u(t) to equation (1.5) with 
5(t ) := S n (t) is a linear combination of v n (t) and w n (t), and more precisely 

u(t) = u'(0)v n (t) + Xu(0)w n (t). 


It follows that 

K(fo )| 2 + A 2 Hf 0 )| 2 < 2|u'(0)| 2 • K(f 0 )| 2 + 2A 2 H0)| 2 • K(t 0 )| 2 

+ 2AV(0)| 2 • |u n (to)| 2 + 2A 4 |u(0)| 2 • \w n (t 0 )\ 2 

= 2|«'(0)p • (K(t»)| 2 + AV(*o)| 2 ) 

+ 2A 2 |m( 0)P • (K(to)p + A 2 |tu„(«o)| 2 ) , 

so that (4.12) follows from (4.13). This proves our original claim. 

If we want a smooth damping coefficient, we need a two steps approximation. First 
of all we choose n 0 € N such that 

I 2 + A 2 k 0 (f 0 )| 2 < e~ 2M , |<(to)| 2 + A 2 |w no (f 0 )| 2 < e~ 2M . 

Then we approximate 8 no (t ) with a nonnegative function 5(t) of class C°° with the 
property that all its time-derivatives vanish both in t — 0 and in t — to- Let v(t) and 
w(t) denote the corresponding solutions to (1.5) with the same initial data of v no (t) and 
w no (t), respectively. From Lemma 4.1 we deduce that, if 5{t) is close enough to 5 no (t) 
in L 2 ((0,fo)), then 

\v'(t 0 )\ 2 + X 2 \v(t 0 )\ 2 < e~ M , \w\t 0 )\ 2 + A 2 |w(t 0 )| 2 < e~ M . 

Since every solution to (1.5) is a linear combination of v(t) and w(t), we can conclude 
exactly as before. □ 

Proof of Theorem 2.1 

Let M be such that 2e -M = e~ Rt °. Let us consider any function 5(t) provided by 
Lemma 4.3 with this choice of M, and let us extend it by periodicity to the half line 
t > 0. Let u(t) be any corresponding solution to (1.5), and let us consider the usual 
energy E(t) := | u’{t)\ 2 + A 2 |u(t)| 2 . 

The estimate of Lemma 4.3 can be applied in all intervals of the form [kto, (k + l)t 0 ], 
yielding that 


E((k + l)t 0 ) < E(kt 0 ) ■ 2e~ M = E(kt 0 ) • e~ Rto Vk e N. 
Therefore, an easy induction proves that 

E(kt 0 ) < E{ 0) • e~ kRt0 Vk e N. 
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Since E(t) is nonincreasing, this implies that 

E(t) < E( 0) exp (— R(t — f 0 ) + ) it > 0, 

which is equivalent to (2.1). 

If we want 5{t) to be piecewise constant, it is enough to take a function with this 
property from Lemma 4.3. If we want 5{t) to be of class C°°, it is enough to take 
from Lemma 4.3 a function of class C°° whose time-derivatives of any order vanish at 
the endpoints of the interval. This condition guarantees that the periodic extension 
remains of class C°°. □ 

Proof of Theorem 2.2 

Let Mk be a sequence of positive real numbers such that 

2e _M ° < ip{2t 0 ), (4.14) 

and 

<p((k + l)t 0 ) • 2e _Mfc < ip((k + 2 )t 0 ) ik > 1. (4-15) 

For every k G N, let 5k : [0, to] y [0,+oo) be one of the functions provided by 
Lemma 4.3, applied with M := M k . Let us define 5 : [0, H-oo) —> [0, +oo) by glueing 
together all these functions, namely by setting 

5(t) := 5 k (t - kt 0 ) it G [kt 0 , [k + l)t 0 ). 

Let u(t) be any corresponding solution to (1.5), and let us consider the usual energy 

E(t) := | u’(t)\ 2 + X 2 \u(t)\ 2 . 

We claim that E(t) < E(0)(p(t) for all t > t 0 , which is equivalent to (2.2). In order 
to prove this result, it is enough to show that 

E(kt 0 ) < E( 0) • (p((k + 1 )t 0 ) ik > 1. (4.16) 

Indeed, since both E(t) and ip(t) are nonincreasing, when t G [kt 0 , (k + l)t 0 \ for some 
k > 1 it turns out that 

E(t) < E(kto) < F'(O) • (p((k + l)t 0 ) < E( 0) • <p(t). 

In order to prove (4.16), we repeatedly apply Lemma 4.3. Since in the interval [0, to] 
the function 5(t) coincides with 5 0 (t), from Lemma 4.3 and (4.14) we deduce that 

E(t 0 ) < E(0) ■ 2e~ Mo < E(0) ■ <p(2t 0 ), 

which proves (4.16) in the case k — 1. 
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Now we proceed by induction. Let us assume that (4.16) holds true for some positive 
integer k. In the interval [kt 0 , ( k + l)t 0 ] the function 5(t) is a time translation of S k (t ), 
and thus we can apply Lemma 4.3 up to this time translation. Keeping (4.15) into 
account, we deduce that 

E((k + l)t 0 ) < E{kt 0 ) ■ 2e~ Mk < E( 0) • <p((k + 1 )t 0 ) • 2e~ Mk < E{ 0) • tp((h + 2 )t 0 ), 
which completes the induction. 

The possibility of choosing a damping coefficient satisfying further requirements 
depends on the analogous possibility in Lemma 4.3 □ 


4.2 Systems of ordinary differential equations 

The following result is the natural generalization of Lemma 4.3 to systems. 

Lemma 4.4 (Decay for all solutions to a system). Let n be a positive integer, let 
(Ai,..., A n ) G (0, Too) n , let M be a positive real number, and let t 0 be defined by (1.8). 

Then there exists a bounded measurable damping coefficient 5 : [0,£ 0 ] [0, +oo) 

such that every solution to (1.7) satisfies 


(K(*o )| 2 + AfcM* 0 )| 2 ) < 

k =1 


E (k(o)l 2 + A> t (o)| 2 ) 

_k= 1 


■ 2e- M . 


(4.17) 


Furthermore, one can choose S(t ) such that 

• either it has the profile shown in Figure 3 on the left, 

• or it is of class C°° and all its time-derivatives vanish both in t = 0 and int = t 0 . 


Proof Let us apply Lemma 4.3 to the k-th equation of the system. Setting t ok := 
7r/(2Afc), we obtain a function 5 k : [0,t 0 fc] —> [0,+cxd) such that every solution to the 
k-th equation of the system satisfies 

K(io*)| 2 + \l\MhO\ 2 < (KMI 2 + Al|« t (0)| 2 ) • 2e- M . 

Now we observe that to — toi + • • • + ton, and we dehne 5 : [0, to] —^> [0, +cxd) by 
glueing together the functions 5 k (t) dehned above. More precisely, we partition [0, to] 
into n subintervals [sfc_i, s k \ with s 0 = 0, s n = t 0 and s k — = t ok , and then we set 


5(t) . d k (t Vt G [sfc—ijS/j). 


Let ... ,u n (t )) be a corresponding solution to the system (1.7), and let us set 


E k (t) := K(t)| 2 + A^Mt)| 2 . 
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We claim that the energy E k (t) of the k-th component is reduced by a factor 2e~ AI 
in the interval [sfc_i,Sfc]. Indeed in this interval 5(t) coincides with S k (t) up to a time 
translation, hence from Lemma 4.3 we deduce that 

Ekito) < Ek(sk) < Ek(sk-i) ■ 2e M < Ek( 0) • 2e M . 

Summing over all indices k from 1 to n we obtain (4.17). 

If we want a damping coefficient with the profile shown in Figure 3 on the left, it is 
enough that all functions Sk{t) are of the form (4.3), with the same value of n for all /c’s 
(this is possible provided that n is large enough). 

If we want a damping coefficient of class C°°, it is enough to choose all functions 
6k{t) of class C°° with all time-derivatives vanishing at the endpoints of the interval. 
This condition guarantees that the glueing procedure yields a function which is still of 
class C°°. □ 

Proof of Theorem 2.3 and Theorem 2.4 

The arguments are analogous to the proofs of Theorem 2.1 and Theorem 2.2, just 
starting from Lemma 4.4 instead of Lemma 4.3. □ 

Remark 4.5. As already mentioned in the introduction, now it should be clear from 
the proofs that the conclusions of Theorem 2.3 and Theorem 2.4 hold true for the 
general equation (1.1), provided that the spectrum of A is finite, even if the dimension 
of eigenspaces is infinite. 

4.3 Partial differential equations 

When the operator A has infinitely many eigenvalues, a constant dissipation is enough to 
dampen all components corresponding to large enough eigenvalues. This is the content 
of next result. 

Lemma 4.6 (PDE with constant dissipation). Let H be a Hilbert space, and let A be 
a self-adjoint nonnegative operator on H with dense domain D(A). Let M be a positive 
real number, and let us assume that 

\A 1/2 x\ 2 > 2M 2 \x\ 2 \/x e D(A 1/2 ). (4.18) 

Then every weak solution u(t) to 

u"{t ) + 2 Mu'(t) + Au(t) = 0 


satisfies 


u'(t )| 2 + \A l/2 u{t)\ 2 < (|m'( 0)I 2 + |A 1/2 w(0)| 2 ) • 8e~ 2Mt Vt > 0. (4.19) 
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Proof Let us consider the usual energy 

m :=i«w + i^y«)i 2 , 

and the modified energy 


E(t):= \u'ft)\ 2 + \A 1 ^ 2 uft)\ 2 + 

2 Mfu! ft), u(t)). 


An easy calculation shows that 



Eft) = -2 ME ft) 

Vt > 0. 

(4.20) 

We claim that 

-| 



-E(t) < E(t) < 2 E(t) 

Vt > 0. 

(4.21) 


Indeed, from the inequality 

2M(u'(t), uft)) < -\u'{t)\ 2 + 2M 2 \u(t)\ 2 
and assumption (4.18) it follows that 

Eft) < Eft) + i|u'(t)| 2 + \A^ 2 u(t)\ 2 < 2 Eft). 

On the other hand, from the inequality 

2 M(uft),u(t)) > ~\uft)\ 2 - |.\/>(/)| 2 
and assumption (4.18) it follows that 

m > E(t) - |i«'(oi 2 - l\A i/2 u(t) r > \E(t). 

From (4.20) and (4.21) we conclude that 

Eft) < 4 Eft) = AE(0)e~ 2Mt < 8E(0)e~ 2M \ 
which is equivalent to (4.19). □ 

The following result is the generalization of Lemma 4.3 and Lemma 4.4 to the infinite 
dimensional setting. 

Lemma 4.7 (Decay for all solutions to a PDE). Let H, A, A& be as in Theorem 2.5. 
Let R be a positive real number, and let 

Tr ' = \ , E f t ( 4 - 22 ) 

A|<2(R+Ai) 2 

Then there exists a bounded measurable damping coefficient 5 : [0, 2 T R ] —> [0, +oo) 
such that every weak solution to (1.1) satisfies 

|u'(2Tr)| 2 + \A 1/2 u( 2T r )\ 2 < (|«'(0)| 2 + \A 1/2 u(0)\ 2 ) ■ e~ R ' 2TR . (4.23) 

Furthermore, one can choose 5ft) such that 

• either it has the profile shown in Figure 4, 

• or it is of class C°° and all its time-derivatives vanish both int = 0 and int = 2 T R . 
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Proof Let us write H as a direct sum 


H — Hr - © H r + , 

where Hr - is the subspace generated by all eigenvectors of A corresponding to eigenval¬ 
ues A 2 < 8(i? + Ai) 2 , and H Rj+ is the closure of the subspace generated by the remaining 
eigenvectors. We point out that Hr - and H R + are ^-invariant subspaces of H. More¬ 
over, the restriction of A to Hr- has only a finite number of eigenvalues, while the 
restriction of A to H R + satisfies the coercivity condition 

\A l/2 x\ 2 > 2(R + Ai) 2 |a:| 2 Vx G D(A 1/2 ) n H R , +t 
and even a stronger condition 

\A 1/2 x\ 2 > 2 (R! + Ai) 2 |a:| 2 Mx G D(A 1/2 ) n H Rt+ , (4.24) 

for some R! > R (actually we can choose 8 (R' + Ai) 2 to be the smallest eigenvalue of A 
greater than 2(R + Ai) 2 ). 

Let u R -(t ) and u R+ (t) denote the components of uft) with respect to the decom¬ 
position, let Eft) := \u'(t)\ 2 + \A 1 / 2 u(t)\ 2 be the usual energy of uft), and let E Rj ±ft) 
denote the energy of the two components. 

Since the restriction of A to H R _ has only a finite number of eigenvalues, the com¬ 
ponent u R — ft) can be regarded as a solution to a system of finitely many ordinary 
differential equations. For this system, the time T R defined by (4.22) coincides with the 
time to defined by (1.8). Therefore, from Lemma 4.4 (see also Remark 4.5) we deduce 
the existence of a function 5 : [0, Tr] —> [0, +oo) which reduces the energy of u R -ft) 
at time t, — Tr by any given factor. In particular, we can choose this factor equal to 
8 e - 2 (i?+Ai)T H an q obtain that 


Er,-(Tr) < Er,_( 0) ■ Se-Wb^. (4.25) 

Now we extend 6 ft) to the interval [0, 2T R ] by setting 6 ft) R + Ai in the second 
half of the interval, namely for every t G (T r ,2Tr\. Since the restriction of A to H R + 
satisfies (4.18), we can apply Lemma 4.6 with M R+ X±. We obtain that 

Er,+ (2T r ) < E Ri+ {T r ) ■ 8e- 2 ( R+Al)T «. (4.26) 

Keeping into account that in both cases the energy is nonincreasing in the whole 
interval, from (4.25) and (4.26) we deduce that 

E(2T r ) < E(0) ■ 8e _2( ' R+Al)Tfi = E(0) ■ 8e~ 2XlTlt ■ e R ' 2TR . (4.27) 

On the other hand, from (4.22) it is clear that 2Ai Tr > n, and hence 8e^ 2Airfl < 1. 
Therefore, now (4.27) reads as 

E(2T r ) < E(0) • e- R2T \ 
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which is exactly (4.23). 

If we want 5(t) with the profile of Figure 4, it is enough to reduce the energy of 
Ur,- (t) through a damping coefficient in [0, Tr] with the profile shown in Figure 3 on 
the left. 

If we want a damping coefficient of class C°° with all time-derivatives vanishing at 
the endpoints, we need an approximation procedure. To this end, we first choose a 
bounded measurable function S(t) in [0, 2 T R ] for which (4.23) holds true for all solutions 
with 2 RTr replaced by a larger constant. This is possible because in the first half of 
the interval we can reduce the energy of UR-(t) by any given factor, and in the second 
half of the interval we can reduce the energy of UR ;+ (t) by a factor e ~ 2( - R+Xi ' )Tr , where 
R' > R is the constant that appears in the reinforced coercivity inequality (4.24) (of 
course we need to set S(t) := R' + Ai in the second half of the period, as required by 
Lemma 4.6). At this point we can approximate 5(t) in L 2 -norm through a damping 
coefficient with the required smoothness, and conclude with the aid of Lemma 4.1. □ 

Proof of Theorem 2.5 

Let us observe that the time 2 Tr, with Tr given by (4.22), coincides with to as defined 
in (1.9). Therefore, from Lemma 4.7 we obtain a function 5 : [0, to] [0, +oo) such 
that every weak solution to (1.1) satisfies 

E(t 0 ) < E(0) ■ e~ Rt °, 

where E(t) := \u'{t)\ 2 + \A 1 ^ 2 u(t)\ 2 denotes the usual energy of the solution. 

Now we extend 8(t) by periodicity to the whole half line t > 0, and we obtain by an 
easy induction that 

E(kt 0 ) < E{0 ) • e~ kRt ° Wk G N. 

Since E(t) is nonincreasing, this implies that 

E(t) < E( 0) • exp (— R(t — f 0 ) + ) Vt > 0, 

which is equivalent to (2.3). 

In analogy with the proof of Theorem 2.1, if we want a piecewise constant damping 
coefficient, we can take a function with this property from Lemma 4.7. If we want a 
damping coefficient of class C°°, it is enough to take from Lemma 4.7 a function of class 
C°° whose time-derivatives of any order vanish at the endpoints of the interval. This 
condition guarantees that the periodic extension remains of class C°°. □ 

Proof of Theorem 2.6 

Let us start by defining (p(t). To this end, let n 0 be the smallest integer such that 
X 2 q > 2X\. For every n > n 0 we consider the positive real number 



and we define T n as in (4.22) with R := R n . Since 2 (R n + Ai) 2 = A 2 , this means that 


T 

1 71. 


n 



l 


Then we set S no -i := 0, C/ no _i := 0, and 


Vn > n 0 . 


n 

Sn := 2 Y T k, 

k=no 


n 

U n := 2 R.T k 

k=no 


for n > no- It is clear that S n and U n are unbounded increasing sequences, with 

2T n+ i = 2 T n + —-< 2 T n + 2 T no < S n 

A n +i 

for every n > Uq + 1, and therefore 


Sn+i — S n + 2T n+ i < 2S n Vn > no + 1. 

Finally, we define (p : [0, +oo) —» (0, +cxd) as the piecewise constant function such 
that 

ip{t) := e~ Un if f 6 [S n , S n+ 1 ) for some n > no — 1. 

It is clear that (pit) is nonincreasing. We claim that it satisfies (2.4). Indeed for 
every n > no — 1 it turns out that 

(p(t)e m = e ~ Un+Rt < e ~ Un+RSn+1 Vt e [S n , S n+1 ), (4.28) 

and for n > no + 1 it turns out that 

n 

—U n + RS n +\ < — U n + 2 RS n = ^ ^ (4i? — 2 R k )T k . 

k=no 

Since R k —> +cxd and T k > T no > 0, we deduce that — U n + RS n+ i —oo. Therefore, 

the right-hand side of (4.28) tends to 0 as n —>■ +oo, which proves (2.4). 

It remains to define S(t). To begin with, for every n > no we apply Lemma 4.7. 
Since the time 2 Tr u of Lemma 4.7 coincides with 2 T n as defined above, from the lemma 
we obtain a damping coefficient 5 n : [0, 2T n ] —> [0, +oo) which reduces the energy of all 
solutions to (1.1) by a factor e~ 2i?7lTn in the interval [0, 2T n ]. 

Now we glue all these functions by setting 

S(t) := S n+ i(t — S n ) if f £ [5 n , S n+ 1 ) for some n > no — 1. 

Let us consider equation (1.1) with this choice of S(t), let u(t) be any solution, and 
let E(t) := \u'(t)\ 2 + \A 1 / 2 u(t)\ 2 be its usual energy. The effect of S(t) in the interval 
[S'n, S'n+i] is the same as the effect of S n+ i(t) in the interval [0, 2T n+1 ], and therefore 

E(S n+ 1 ) < E(S n ) • e ~ 2Rn + lTn + 1 Vn > n 0 - 1. 
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At this point an easy induction shows that 

E(S n ) < E( 0) • exp (-2 (R no T no + ... + R n T n )) = E(0) ■ e~ Un Vn > n 0 . 

On the other hand, also for n = no — 1 it is true that E(S n ) = E( 0) • e~ Un , for 
the trivial reason that f7 no _i = 0. Since E(t) is nonincreasing for i > 0, we can finally 
conclude that 

E(t) < E(S n ) < E(0)e~ Un = E( 0) • f{t) 

for every n > no — 1 and every i G [S n , S^+i], which proves (2.5). □ 

4.4 Further requirements on the damping coefficient 

Proof of Proposition 2.7 

Let us consider the usual energy E(t) := \u'(t)\ 2 + A 2 |u(i)| 2 . Its time-derivative is 
E'(t) = —45(t)\u'(t)\ 2 , so that 

E\t) > -4 5(t)E(t). 

Integrating this differential inequality we deduce that 

E(t) > E(0) ■ exp ^—4 J S(s) ds^j Vi > 0. 

This is enough to deal with the first two cases. 

For the third statement, let us assume now that 5(t) > A for every t > T. In this 
case we consider the Riccati equation 


<p'(t) = A 2 - 2 6(t)(p(t) + 


(4.29) 


Let us set T* := max{T, 1/A}. It is not difficult to check that (p(t) = 0 is a subsolution 
of (4.29) for t > 0, while (p(t) := A — 1/i is a supersolution for t > T* owing to our 
assumptions on S(t). As a consequence, the solution to (4.29) with </?(T*) = 0 is dehned 
for every t > T* and satisfies 


0 < (f{t) < A — - Vi > T* 


Finally, one can check that 

u(t) := exp 


f(s) ds 


Vi > T* 


(4.30) 


defines a solution to (1.5), which can be easily extended to the whole half line t > 0. 
Due to the estimate from above in (4.30), this solution satisfies 


u(t)\ > cte xt Vi > T* 
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for a suitable positive constant c. Since equation (1.5) is linear, c 1 u(t ) is again a 
solution, and satisfies (2.6). □ 

The rest of this section is devoted to the proof of Theorem 2.8. This requires two 
preliminary results. In the first one we show how to dampen one solution to (1.5) 
through a slowly increasing damping coefficient. 

Lemma 4.8 (Supercritical energy reduction of a solution). Let A and M be positive real 
numbers, and let e e (0, A). Let us define 


M 2 

h: ~ fe + V 

S(t):—X + £t Vt e [0,ti]. 

Then there exists a nonzero solution v(t) to (1.5) such that 

|u'(t 1 )| 2 + A>(f!)| 2 < (|u'(0)| 2 + A>(0)| 2 ) e~ Mtl . 

Proof Let us start by proving that (4.31) implies that 

(3A 2 + 8e 2 t 2 ) exp(—2Ati — 2 et\) < 2A 2 exp(— Mtf). 
Indeed from (4.31) it turns out that 


2 £t{ > Mt i + —t i. 

A 


2Ati > 4, 2 £t{ > Mt i, 

From the first two inequalities it follows that 

3A 2 exp(—2Ati - 2 et\) < A 2 • 3e" 4 • e~ Mtl < X 2 e~ Mtl 
From the third inequality it follows that 

exp(-2Ati - 2 etl) < e" 2£ ‘i < e “ 4etl/A • e~ Mtl . 
Therefore, from the inequality 


Ae 


t 2 e~ ct < 


e 2 c 2 


Vt >0, Vc > 0, 


applied with c := Ae/X and t := f l5 we obtain that 

8 £ 2 t\ exp(—2Ati — 2 et\) < %e 2 t\ exp(—^ A£t\/X) ■ exp (—Mt\) 


< 8e" 


4A 2 


16e 2 £ 2 


exp (-Mt! 


< A 2 exp {-Mti). 


(4.31) 

(4.32) 

(4.33) 


(4.34) 


(4.35) 


(4.36) 
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Adding (4.35) and (4.36) we obtain exactly (4.34). 

We are now ready to prove the existence of v(t). Let us consider again the Riccati 
equation (4.29). When S(t ) is given by (4.32), a simple computation shows that (p(t) := 
A + 2 et is a supersolution for t > 0, and the constant function ip{t) := 2(A + 2et\) is a 
subsolution for t G [0,£i]. ft follows that the solution to (4.29) with “final” condition 


— A + 2sti 

is defined for every t G [0, 1 1 ] and satisfies 

(pit) > A + 2 et Vf G [0, t\]. 
At this point a simple computation shows that 


v(t) := exp ( — / ip(s) ds 

i o 


V£ G [0, t\ 


is a solution to (1.5). We claim that this solution satisfies (4.33). 

To begin with, we observe that v(0) = 1, and from (4.38) it follows that 

v(ti) = exp f ip(s) ds^j < exp (—Ati — etf) . 

As for the time-derivative, it is given by 


v'(t) = —(f(t) exp ( — / (p[s) ds ) , 


hence by (4.38) 

and by (4.37) and (4.38) 


K(0)| = M0)|>A, 


|u ; (£i) | = \(f{ti) | exp J ip(s) ds^j < (A + 2et\) exp (—A£i — etf) . 

Recalling (4.34), from all these estimates it follows that 

|r/(£i)| 2 +A 2 |u(£i)| 2 < [(A + 2e£i) 2 + A 2 ] exp (—2A£i — 2e£ 2 ) 

< [3A 2 + 8 e 2 tf) exp (—2Afi — 2 etf) 

< 2A 2 exp(— Mti) 

< (|^'(0)| 2 + A 2 |u(0)| 2 )exp(-M£ 1 ), 


(4.37) 


(4.38) 


which proves (4.33). □ 

The second preliminary result clarifies that any constant damping coefficient below 
the threshold A allows solutions to rotate in the phase space. 
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Lemma 4.9 (Subcritical rotation in the phase space). Let X be a positive real number, 
and let e 6 (0, A). Let (a, j3) and ( 7 , 5) be two pairs of real numbers with a 2 + /3 2 7 ^ 0 
and 7 2 + 8 2 7 ^ 0 . 

Then there exist positive real numbers r and t 2 , with 

2?r 

t 2 < , 

yj e( 2A — e) 

with the following property. The solution u(t) to equation (1.5) with constant dissipation 
8(f) := A — e and initial data w(0) = a and r/(0) = f3 satisfies 

u(t 2 ) = T 7 , u'(t 2 ) = rS. 


(4.39) 


Proof The solution u(t) can be written in the form u(t) := e~^ £ ' t v(t), where v(t) is 
a solution to 

v"(t) + e(2X — e)v(t) = 0. 

Integrating this differential equation we see that the pair (v(t),v'(t)) rotates in the 
phase space with period equal to the right-hand side of (4.39). During the complete 
rotation the pair (v(t),v'(t)) turns out to be a positive multiple of any given vector. 
Since the same is true for u(t), this proves the result. □ 


Proof of Theorem 2.8 

Strategy Let us describe the strategy of the proof before entering into details. Let 
us consider the function 8(t) with the graph as in Figure 5, and then extended by 
periodicity. 



Figure 5: profile of 8(t) for the proof of Theorem 2.8 

The value at the endpoints is A, the value in the horizontal part is A — e, the slope in 
the oblique sections is ±e. This graph depends on two parameters to be fixed, namely 
the length t\ of the first and last oblique sections, and the length t 2 of the horizontal 
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plateau. The length of each of the two oblique sections attached to the horizontal part 
is necessarily + 1. The total length is therefore 


to :=4f 1 + f 2 + 2. (4.40) 

Now the idea is the following. In the first ascending section, 5(t) causes a great 
reduction of the energy of a special solution v(t) to (1.5). Let w(t) be the solution with 
initial data orthogonal to the initial data of v(t). In the first two oblique sections we 
know that the energy of w(t) is just nonincreasing. On the other hand, we can choose 
t ,2 in such a way that, at the end of the horizontal plateau, w(t) has the right “initial” 
data which guarantee a reduction of its energy in the third oblique section. 

In other words, the first ascending section cuts the energy of v(t), while the horizontal 
plateau rotates w(t) preparing it to undergo a cut of its energy during the second as¬ 
cending section. The two descending sections are merely junctions between the “active” 
parts of the graph. 

Choice of parameters Let us set 

M := 8(7t + e)R + 4e log 2 + 1, (4.41) 

and let ti be defined by (4.31). The function S(t) we are considering is equal to A + et 
for every t G [0,£i]. Therefore, from Lemma 4.8 we know that (1.5) admits a solution 
v(t) satisfying (4.33). Since the equation is linear, we can always assume that 

|i/(0)| 2 + A>(0)| 2 = 1. (4.42) 

In the sequel we also need to consider 8(t) as extended to t G [—1, 0] with the same 
expression A + et. Accordingly, we can extend v(t) to the interval [—1, 0], and consider 
the pair (v(—1), v'(— 1)) of its data for t = —1. 

Let us consider now the solution w(t) to (1.5) with initial data 

w(0) = w'(0) = —Au(0). 

A 

We point out that the initial data of w(t) satisfy 

|w'(0)| 2 +A 2 |w(0)| 2 = 1 (4.43) 

and they are orthogonal to the initial data of v(t) in the sense that 

w'(0)u'(0) + A 2 u?(0)u(0) = 0. (4.44) 

We are now ready to choose t 2 - Let us consider the end of the second oblique section, 
corresponding to t — 2 1\ + 1, and let us set 


(a,/3) := (w(2t 1 + l),w'(2t 1 + 1)), 


(7,5) := (u(-l),t/(-l)). 



From Lemma 4.9 we deduce the existence of a positive time t 2 satisfying (4.39) such 
that the effect of the constant dissipation equal to A — £ on the solution with “initial” 
data (ct,/3) is to transform it in a multiple of (7,5) in a time t 2 . It follows that, at the 
end of the horizontal plateau at time t — 2ti + t 2 + 1, the solution w(t) is in the same 
conditions as (the extension of) the solution v(t) at time t — — 1. As a consequence, 
the energy of w(t ) is reduced by a factor e~ Mtl during the third oblique section, from 
t = 2ti + t 2 + 1 to t = 3ti + t 2 + 2. 


Estimates From (4.40), (4.31) and (4.39) it follows that 


2 M 8 „ 

to — -b -r + 2 + t 2 < 

£ A 


2 M 8 „ 2tt 
-+ +2 + —, 

£ A £ 


so that (2.7) follows from (4.41). 

As for energy estimates, let E u (t ) := \u'(t)\ 2 + X 2 \u(t)\ 2 denote as usual the energy of 
a solution to (1.5). The energy of v(t) has been reduced in the first ascending section. 
Since it is always nonincreasing, it follows that 


E v (t 0 ) < E v (t i) < E v { 0) • e~ Mtl = e~ Mtl . (4.45) 


Similarly, the energy of w(t) has been reduced in the second ascending section, hence 
E w (t 0 ) < E w (3ti + t 2 + 2) < E w (2ti + t 2 + 1 ) • e~ Mtl < E w ( 0) • e~ Mtl = e~ Mt \ (4.46) 


Let us consider now any solution u(t) to (1.5). Since v(t) and w(t) are linearly 
independent, we can write 

u{t) = av(t) + bw{t) 

for suitable real constants a and b. Thanks to the orthonormality relations (4.42), (4.43), 
and (4.44), it turns out that 

a 2 + b 2 = E u ( 0). 

Moreover, from (4.45) and (4.46) it follows that 

E u (t 0 ) = | an'(t 0 ) + bw'(t 0 ) | 2 + A 2 |au(t 0 ) + bw(t 0 ) | 2 

< 2a 2 E v (to) + 2b 2 E w {to) 

< 2 (a 2 + b 2 )e~ Mtl 
= E u (0) ■ 2e~ Mtl . 

Now we claim that 

2e~ Mtl < e~ Rt ° , (4.47) 

and hence 

E u (t 0 ) < E u ( 0) • e' ra °. 

Indeed, from (4.41) it follows that M > 8 R, and hence 

^Mti > 4Rt u (4.48) 
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while from (4.41) and (4.31) it follows that 


1 M 2 M ( 2vr\ „ , 

-Mti > + — J R + log 2. 


2 4e 4e 

Summing (4.48) and (4.49) we obtain that 

2t r\ 


(4.49) 


Mt\ > ( 4 ti + 2 H-J R + log 2 > (4 1\ + 2 + t^R + log 2 — 4“ log 2, 


which is equivalent to (4.47). 

Finally, when we extend 5(t) by periodicity to the whole half line t > 0, an easy 
induction gives that 

E u (nt 0 ) < E u { 0) • e~ nRt0 Vn e N. 

Since E u (t) is nonincreasing, this implies (2.8). □ 


5 Applications to PDEs 

In this section we present some simple examples of application of the abstract results 
stated in section 2.3. 

Let 11 C M. d be a connected bounded open set with smooth boundary. We consider 
two model examples: a dissipative wave equation 


u tt — A u + p(x)u + 5(t)u t = 0 

t > o, iefl, 

(5.1) 

and a dissipative beam/plate equation 



u tt + A 2 w + q{x)u + 8(t)u t = 0 

t > 0, x G fl, 

(5.2) 


where p(x) and q(x) are nonnegative bounded measurable functions. Just to fix ideas, 
we consider equation (5.1) with homogeneous Dirichlet boundary conditions u — 0, and 
equation (5.2) with one of the following boundary conditions: either u = A u = 0 (simply 
supported beam or plate) or u — \ Vu | = 0 (clamped beam or plate). 

In both cases Theorem 2.5 and Theorem 2.6 imply that 

• any exponential decay rate e~ Rt can be realized through a suitable periodic damp¬ 
ing coefficient with period Tr that depends on R , 

• with a suitable non-periodic coefficient one can achieve a decay rate (pit) faster 
than all exponentials. 

The period T R and the ultra-exponential decay rate ip{t) depend on the growth of 
eigenvalues, which in turn is known to depend on the space dimension d. More precisely, 
when eigenvalues are arranged in increasing order, the n-th eigenvalue is comparable 
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with the power of n indicated in the hfth column of the table below (for a proof we refer 
to the seminal papers [1, 7, 10]). 

In the third and fourth column we state our estimates for Tr and <p(t). They have to 
be interpreted as asymptotic behaviors, namely when we write log-R in the third column 
we actually mean that Tr = O(log-R) in that case. In the same way, we write shortly 
exp(— t a ) instead of exp(— 0(t a )). 


Equation 

d 

Tr 

(p(t) 


T 

± n 

Sn 

u n 

(5.1) 

1 

log R 

exp(— t 2 / logf) 

n 

log n 

n log n 

n 2 log n 

(5.1) 

> 2 


exp(— t 2d T 2d ~ 1 '>) 

n l / d 

n l ~ l/d 

n 2 ~ l/d 

n 2 

(5.2) 

1 

1 

any 

n 2 

— 

— 

— 

(5.2) 

2 

log-R 

exp(— t 2 / log t) 

n 

log n 

n log n 

n 2 log n 

(5.2) 

> 3 

Rd/2-l 

exp (—tVG-L) 

n 2 / d 

n l ~ 2/d 

n 2 ~ 2/d 

n 2 


In particular, we observe that for equation (5.2) with d — 1 (beam equation), the 
series of reciprocals of eigenvalues is convergent. As a consequence, the period Tr given 
by (1.9) is bounded independently of R. At this point, the same argument of the proof 
of Theorems 2.2 and 2.4 gives that any decay rate <p(t) can be achieved if we are allowed 
to exploit non-periodic damping coefficients. 

In all other cases, the series of reciprocals of eigenvalues is divergent, and Tr grows 
with R according to (1.9). In order to compute (p(t), we need to estimate the growth 
of the sequences S n , T n , U n defined in the proof of Theorem 2.6, and recall that ip(t) 
decays in such a way that (p(S n ) = exp(— U n ). The computations are reported in the 
last four columns of the table, with the usual agreement that entries have to be intended 
in the sense of the “big O” notation. 

The optimality of the third and fourth column of the table might probably deserve 
further future investigation. 
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